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Abstract. We study the luminosity function of the high-redshift galaxy population with redshifts 3 < z < 4 using a purely 
/-band magnitude- selected spectroscopic sample obtained in the framework of the VIMOS VLT Deep Survey (VVDS). We 
determine the luminosity function from the VVDS, taking care to add as few assumptions and as simple corrections as possi- 
ble, and compare our results with those obtained from photometric studies, based on Lyman-break selections or photometric- 
redshift measurements. We find that in the redshift range 3 < z < 4, the VVDS luminosity function is parameterized by 
</>' = 1.24^Q"jg 10"^ mag"' Mpc"^ and M* = -21.49^Q"[g, assuming a slope a = -1.4 consistent with most previous studies. 
While (p' is comparable to previously found values, M' is significantly brighter by about 0.5 mag at least. Using the conservative 
slope a = -1.4, we find a luminosity density at 1700 A £nao(M < -18.5) = 2.4 10" W Mpc^^ and £]°'^l, = 3.1 10" W Mpc^^, 
comparable to that estimated in other studies. The unexpectedly large number of very bright galaxies found in the VVDS indi- 
cates that the color-selection and photometric-redshift techniques that are generally used to build high-redshift galaxy samples 
may be affected by a significant fraction of color-measurement failures or by incomplete modelling of the mix of stellar emis- 
sion, AGN contribution, dust absorption and intergalactic extinction assumed to identify high-redshift galaxies, making pure 
magnitude selection better able to trace the full population. Because of the difficulty to identify all low-luminosity galaxies in 
a spectroscopic survey, the luminosity density could still be significantly underestimated. We also find that the relative contri- 
bution of the most luminous galaxies compared to the fainter ones is at least twice as large in the VVDS compared to former 
estimates. Therefore, the VVDS paints a quite different picture of the role of the most actively star-forming galaxies in the 
history of star formation. 
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1. Introduction 

Galaxies are formed and evolve through the influence of com- 
plex physical processes. Galaxy number counts are one of 
the most obvious observational signatures of these processes, 
and constitute, as a consequence, a key ingredient to their un- 
derstanding. Notwithstanding how galaxies are assembled, the 
more galaxies there are, the more efficient star formation at ear- 
lier epochs must have been. Counting the numbers of galaxies 
as a function of their total stellar luminosity and as a function 
of redshift in an exhaustive manner is therefore a necessary pre- 
requisite to the determination of the total number of stars at a 
given epoch in the life of the universe. 

Star-forming galaxies are most conspicuously revealed in 
the ultraviolet domain. The building of luminosity functions 
(LFs) in the ultraviolet, which describe the galaxy number den- 
sities as a function of their ultraviolet luminosities, is there- 
fore an important step in the study of the star-formation his- 
tory. At redshifts above two the detection of these galax- 
ies is in theory made simpler by the fact that ultraviolet 
photons are redshifted in the optical range, which can be 
easily observed. However, building large samples of high- 
redshift galaxies is challenging. At redshifts z » 1 the huge 
amount of observing time required to measure large num- 
bers of redshifts has forced up to now to apply preselection 
criteria to weed out the much more numerous low-redshift 
galaxies. A number of deep redshift surveys have identi- 
fied various types of galaxies at redshifts above 2, includ- 
ing Lyman-b reak galaxies (LBGs; Steidel et al. 19961) 



distant 
axies 



red galaxies (van Dokkum et al. 20031) . BzK-selected gal 
JPad di et al 
1200^ 



2004?), faint K-selected galaxies (Abraham et al.l 
2004), or counterparts o f faint sources detect ed in radio, X- 
ray or sub-mm surveys dChapman et al.ll2005b . In addition, a 
number of photometric-redshift surveys have identified high- 
redshift galaxies based on the expected spect rophotometric 
signature of combined stellar populations (e.g., 'pontana et al.' 
2000; Polietal. 2001; Arnouts et al. 2002; Rowan-Robinson 
20031; iGabasch et al.1 l2004t iSaracco et all 12006 ). The under- 



standing of the selection eff'ects is however a serious challenge 
for the making of a complete census of the high-redshift galaxy 
population. In particular, the different criteria used to pick-up 
high-redshift galaxies ha ve been shown to p roduce distinct, but 
overlapping populations (iReddv et al ]l2005h . 

In spite of these difliculties, several groups have buil t ultra- 



violet luminosity functions at redshift 3 < z < 4. Steide l et al 
d 1 9991) )uilt the first ultraviolet LFs at z ~ 3 and z ~ 4 based on 



the first large sample of high-redshift galaxies ph otometrically 
selected using the LBG color-color d iagrams ( Steidel et al.l 
I996h. Using a very similar approach, Sawicki & T hompson 



( l2006ah extended the study of the ultraviolet LF of LBGs using 



* Based on observations collected at the European Southern 
Observatory Very Large Telescope, Paranal, Chile, program 070.A- 
9007(A) 

** Based on observations obtained at the Canada-France-Hawaii 
Telescope (CFHT), which is operated by the National Research 
Council of Canada, the Institut National des Sciences de I'Univers 
of the Centre National de la Recherche Scientifique of France, and the 
University of Hawaii. 



the small, b u t very deep. Keck Deep Field (KDF). At z ~ 4, 



Ouchi et al.l (|2004|) adapted the LBG selection technique to 
the filters used in the Subaru Deep Survey (SDS). Luminosity 
functions based on photometric-redshift studies have been cal- 
culated from samples built usi ng deep imaging surveys, like 



die H ubble Deep Field (HDF; iPoU et a l. 2001; Arnouts et al 



20051) and the Fors Deep Field (FDF; IGabasch et al..,2004t) . It 



must be pointed out that the difficulty in obtaining large sam- 
ples without preselection have prevented up to now the study 
of high-redshift ultraviolet LFs based on spectroscopic sam- 
ples, although spectrosc opic confirmations have been obtained 



for the LBG sample of iSteidel et al.l (Il996l) 



The VIMOS VLT Deep Survey (VVDS; i Le Fevre et al 



2005bl) is following a unique approach combining both a 



deep magnitude-selected spectroscopic survey and a large sam- 
ple, an approach made possible by the recent appearance 
of the VIMOS high-multiplex multi-object spectrographs on 
the European South ern Observatory Very Large Telescope 
dLe Fevre et al.ll2003l) . The /-band selection of objects brighter 
than magnitude Iab - 24 allows us to identify galaxies up to 
z - 5, and in the "First Epoch" VVDS we find 970 candidat e 
galaxies with a redshift 1.4 < z < 5 (Le Fevre et al. 2005al) . 



A surprising result is that the pure magnitude selection identi- 
fies a n unexpectedly large p opulation of galaxies at high red- 
shift ( Le Fevre et al ] |2005al) . with surface densities of galaxies 
at redshifts z ~ 2 to z ~ 4 two to six times higher than reported 
by previous studies using LBG selection techniques. 

In this Paper, we aim to quantify the luminosity function of 
the population of galaxies at high redshift found in the VVDS, 
concentrating on the 3 < z < 4 redshift range. We study the 
implications of the ultraviolet luminosity function on the ultra- 
violet luminosity density at these redshifts. The evolution of the 
ultraviolet lumi nosity density and t he history of star formation 
are presented in Tresse et al. ( 20061) . 

Throughout the paper we use the concordance cosmology 
with //() = 70 km s"' Mpc-^ Qm = 0.3 and Qa = 0.7. 



2. Data 

We use here the VVDS "Fi rst Epoch" sample in th e VVDS-02h 
(-H02h 26m -04° 30') field dLe Fevre et al.l2005bl) . This sample 



consists of the spectra of 9295 objects obtained over a sky area 
of 1720 arcmin^ using the VIMOS multi-obje ct spectrograph 



on the VLT unit 3 Melipal at ESO in Paranal dLe Fevre et al 
E003) . Targets between limiting magnitudes /ab - 17.5 and 
/ab = 24 have been randomly selected. On average, spectra 
have been obtained for 24% of the sources in the photometric 
catalogue in the above magnitude range. The VIMOS observa- 
tions have been carried out using the low-resolution red grism, 
which provides spectra with a resolution ^ 230 over the 5500- 
9500 A wavelength range. Typically 500-600 spectra are ob- 
tained simultaneously. The nominal exposure time is 16200 s, 
split into 10 separate exposures of 27min, shifting the objects 
along the slits to reduce the effect of fringing at long wave- 
lengths. Data reduct ion has been perform ed with the VIPGI 
data processing tool ("Scodegg io et al.l2005 ). Further details are 
given in .Le Fevre et al., (,2005blK 
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As explained in lLe Fevre et al. I (l2005bh . each spectrum has 
been assigned a flag which indicates the reliability of the red- 
shift determination. Objects with flags 3 and 4 have very re- 
liable measurements. A flag 2 indicates a reasonably reliable 
measurement, but with a non-negligible probability that the 
identified redshift is wrong; a flag 1 indicates a tentative red- 
shift, with a significant probability of misidentification. When 
no redshift could be meaningfully estimated from a spectrum, 
the object get s a flag and no redshift determination. In 



Le Fevre et alJ (I2005M) . confidence levels of 50%, 80%, 95% 
and ~ 100% have been estimated for the flag 1, 2, 3 and 4, re- 
spectively, on the basis of repeated observations; these levels 
have ho wever been determined mostly for low-redshift galax- 
ies, and lLe Fevre et al.l (12005 ah refined these estimates for spe- 
cific redshift ranges at z > 1 .4. In total, among the 266 objects 
which have been attributed a redshift between 3 and 4, 12 ob- 
jects have been attributed a flag 3 or 4, 101 a flag 2 and 153 
a flag 1. In addition, 665 objects have a flag 0, and some of 
them may possibly be located at redshifts 3 < z < 4. Objects 
classified as QSOs are excluded from our sample. 

To build the luminosity functions, we complement our 
spectroscopic sample with multi-band photometry. The VVDS 
spectroscopic sample has been set up using photometric ob- 
servations in the B, V, R and / bands obtained with the 
CFH12K CCD camera on the Canada-France-Hawaii tele- 



scope (CFHT) dMcCracken et alj|2003l: iLe Fevre et alj|2004l) . 
In order to extend the photometric coverage down to the U 
band, which is very important for the study of galaxies at 
z > 3, we use here deep u*, g', r', i' and z' photometry 
obtained in the framework of the Legacy Survey project at 
the CFHT (CFHTLsQ) with the MegaPrime/MegaCam imager. 
Processing of the CFHTLS photometry is performed by the 
TERAPIX consortiunfl CFHTLS's Dl deep field covers the 
VVDS-02h deep field. We use here the T0002 release, whose 
accumulated exposure times are, for the u*, g',r', i' and z', 
9.9hr, 5.4hr, 14.6hr, 33 hr and 15hr respectively. The follow- 
ing 50%-completeness limiting magnitudes have been obtained 
in the individual bands: m*,,^^ = 26.4, g'^^^^ - 26.3, r^g^^ - 26.1, 
'50% ~ ^^-^ ^^'^ 4o% ~ 24.9. The CFHTLS photometric cat- 
alogue has been constructed using the VVDS astrometric ref- 
erence fra me, which shows excellen t relative astrometric per- 
formances dMcCracken et al.l l2003h . Therefore only a small 



matching aperture of ^ 0.3 arcsec was needed when matching 
the VVDS spectroscopic catalog with the CFHTLS photomet- 
ric catalog, avoiding multiple identifications. 

3. Weights and corrections in the calculation of 
the luminosity function 

While the VVDS selection function is very simple due to the 
VVDS being a purely magnitude-limited survey, the fraction of 
unobserved objects and errors in redshift determinations must 
be taken into account in the calculation of the luminosity func- 
tions. A specific weight is therefore applied to each object in 
order to recover the proper volume density at a given luminos- 



ity. The weight includes three effects: the probability of a given 
object to have been observed spectroscopically; the probability 
of a given redshift to be correct; and the fraction of objects with 
a redshift in the range 3 < z < 4 mistakenly classified at other 
redshifts. 

3.1. Comparison with previous LF studies witiiin tine 
VVDS 

These difficulties have already been addressed in the frame- 
work of the VVDS for the determination of LFs at redshifts 



z < 2 by lllbert et al.l (120051), and the rnethodo lo gy has been 
reused in subsequent works (lllbert et al. I l2006bll5 IZucca et al 



.2006) . Except for the effect of probability of observation we 
use here a different approach from that used in these papers. 



In lllbert et al.l (l2005h . flag-2, flag-3 and flag-4 objects have 
been assumed to have correct redshift determinations, while 
flag-1 objects have not been used. As already mentioned, the 
fraction of correct redshifts for flag-2 objects is about 80% at 
low redshift. The ~ 20% of wrong redshifts have however neg- 
ligible impact on the determination of the LF, since, at low 
redshift, flag-2 objects constitute only about 40% of the ob- 
jects with good-to-excellent redshift determination. At redshift 
3 < z < 4 this figure reaches 90%. Moreover, at these red- 
shifts, the fraction of correct redshifts among the flag-2 ob- 
i ects is markedly lowe r than 80%, and actually closer to 50% 
( Le Fevre et al. I l2005ah . As a consequence, we cannot assume 
that all flag-2 objects have correct redshift determination, since 
the wrong-redshift objects would be almost as numerous as 
those with correct redshifts. 



The discarding of all flag-1 objects in lllbert et al.l (120051) 
makes that they must be tre ated as an i ncomp leteness, sim- 



ilarly to the flag-0 objects, lllbert et al.l (l2005l) summed the 



photometric-redshift probability densities for all flag-0 and 
flag-1 objects to determine the redshift distribution of these ob- 
jects, which allowed them to estimate the fraction of missed 
objects (i.e., having received either a flag or a flag 1) as a 
function of redshift and magnitude. Again, here we have to opt 
for a different approach, because photometric redshifts are not 
reliable enough at redshift z > 2 to allow a proper estimate of 
the correction to be applied. 

Below we explain in detail our treatment of the different 
effects. 

3.2. Sampling correction 

The target sampling rate (TSR) is the probability that a given 
object in our field has been observed spectroscopically with 
VIMOS. The tool to pla ce the shts for flie VVDS Deep Field 



' http ://w w w. cfht . hawaii .edu/S cience/CFHLS/ 
^ http://terapix.iap.fr 



observations (VMMPS; ' Bottini et al.l 12005!) has been used in 
"optimized mode", which provides the largest multiplexing 
possible, instead of the "random mode", where all objects have 
the same probability to be observed. This has the consequence 
that the TSR of a given source depends on the size of the ob- 
ject, because extended objects have so mewhat smaUer cha nce 
of being targeted than point-like objects, lllbert et al.l (120051) es- 
timated the TSR as a function of the X-radius parameter, which 
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is the size of the object along the slit axis. We use here th e same 
estimates, which are reported in figure 1 of lllbert et al. 1 12005). 
To take this effect into account, the weight of source / is multi- 
plied by 1 ITSRi. On average, the TSR is of the order of 24%, 
and the effect of the X-radius parameter is small, the sources 
being mostly point-like. 



3.3. Wrongly assigned redshifts 

The galaxy counts must be corrected to account for the frac- 
tion of objects with w rong redshift determination. Following 



Le Fevre et al.l (l2005ah . we try to estimate for each flag cat- 
egory the correct-redshift fraction (CRF), i.e. the fraction of 
objects with correct redshift determination, by building the 
stacked composite spectrum of all flag-3/4 objects (whose 
CRFs are assumed to be equal to 1), which we compare with 
stacked composite spectra of flag-1 and flag-2 objects respec- 
tively. The equivalent widths of the different spectral lines are 
then compared to those found in the reference spectrum con- 
sisting of only flag-3/4 objects. Assuming that the population 
of galaxies is similar for each flag, the average ratio between 
the equivalent widths in the less reliable-redshift stacked spec- 
tra and the reference spectrum provides an estimate of the CRF 
for the objects with quality flags 1 and 2. We find a CRF of 
0.44 for the flag-1 objects and of 0.49 for the flag -2 objects, 
consistent with the results of iLe Fevre et al.l (l2005al) . and again 
much lower than the value of 80% found at lower redshift 
jLe Fevre et al.ll2005bh . 

There are two ways to use the CRFs. The simplest assump- 
tion is that all sources in a flag category have the same proba- 
biUty to have a correct redshift. This corresponds to multiply- 
ing the weight of a given source with a flag k by CRFk, CRFk 
being equal to 0.44 for k - \ and to 0.49 for k - 2; we re- 
fer to this method as the equal weight. Alternatively, as a red- 
shift is either correct or wrong, we can try to identify objects 
that have smaller chances of having correct redshifts and dis- 
card them. The " Algorithms for Luminosity Function" (ALF; 



Ilbert et al.ll2005h . which is the tool developed in the frame- 



work of the VVDS consortium to calculate LFs, provides us 
with a quantitative criterion: When calculating absolute lumi- 
nosities, ALF determines a merit figure based on the agree- 
ment, at the given redshift, between the best fitted template and 
the 9 photometric filters used here. We can therefore, for flag-1 
and flag-2 categories, select the CRFi^ x {k = 1,2) objects 
having the best merit figures, where Nt is the number of objects 
with flag k, and discard the remaining objects. We refer to this 
method as the photometric rejection. We point out that mixing 
separate photometric sets can lead to problems with discrepant 
0-points (in fact, th i s can h appen even within a photometric 
svstem). [ilbert et aL (2006a) determined that the most accurate 
photometric redshifts that can be obtained combining CFHTLS 
and CFH12K photometry requires adjustement in the 0-points 
of the different filters. However, these adjustm ents are always 



small er than 0.1 mag for the CWW templates (IColeman et al 



1983); we therefore neglect this correction here. 




1200 



1400 



1600 



1800 



Wavelength 



Fig. 1. Composite spectra of selected candidate 3 < z < 4 ob- 
jects. Top panel: Composite spectrum of the 12 flag-3 and flag- 
4 objects. Middle panel: Composite spectrum of the 117 flag- 
1 and flag-2 spectra selected using the photometric-rejection 
method; flag-3 and flag-4 objects are not included. Bottom 
panel: Composite spectrum of the 137 flag-1 and flag-2 spectra 
rejected by the photometric-rejection method. The major lines 
are indicated. The normalization has been chosen so that the 
equivalent widths of C iv can be visually compared. 



Fig. [U shows the composite spectra of the 12 flag-3/4 ob- 
jects, the 117 flag-1/2 selected spectra and of the 137 flag-1/2 



rejected ones. The composite spectrum of the selected objects 
shows all the expected features of galaxies at z ~ 3-4, both 
in the spectral lines and in the continuum. While the com- 
posite of rejected spectra stills shows some of the lines ex- 
pected in high-redshift galaxies (especially Civ /1 1549 A and 
Si 11/11527 A), they are significantly weaker and some strong 
lines, like Feii/11609A, are missing. Part of the lines appear- 
ing in this composite might therefore be the result of wrong 
identification of features in the spectra with Civ or Sin and 
possibly another shorter-wavelength line. In addition, the con- 
tinuum shape is not really as expected, especially around Ly a. 
We conclude that the photometric-rejection method very signif- 
icantly increases the fraction of correctly-identified galaxies in 
our sample, and while it may be that some correctly-identified 
galaxies are improperly rejected, the remaining contamination 
is probably quite small. 
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3.4. High-redshift galaxy incompleteness 

Incompleteness affects all types of surveys. For instance, to de- 
termine U V luminosity function at high redshift through LBG 



selection dSteidel et alj 119991: ISawicki & ThompsonI l2006al 3 



hereafter S99 and ST06A respectively), the fraction of galax- 
ies escaping the selection needs to be estimated using com- 
plex simulations of the effects of photometric incompleteness, 
unadapted flux integration aperture, contaminations, statisti- 
cal scatter in the measured magnitudes, etc. Thanks to the 
very simple selection function of the VVDS, such simulations, 
which rely on assumptions that are not always easy to assess, 
are not necessary. However, incompleteness does arise from 
the redshift-determination process through the lack of sufficient 
spectral information or the misidentification of the spectral fea- 
tures. 

The redshift of some objects, especially the faintest ones, 
could not be determined; these objects are the so-called flag-0 
objects, and contribute directly to the VVDS incompleteness. 
Assuming that the redshift distribution of these objects is not 
very different from that of the rest of the population, we can 
correct our galaxy counts to take these flag-0 objects into ac- 
count: We estimate the fraction of flag-0 objects as a function of 
/ab magnitude. This fraction, /o(/ab), increases steeply at faint 
magnitudes. Thus, for each galaxy / in the range 3 < z < 4 with 
a magnitude /ab,, we apply a correction factor 1/(1 -/o(/ab,))- 
This correction will have a negligible impact on the bright part 
of the LF, because the fraction of flag-0 objects is very low at 
bright magnitudes . 

In 'l lbert et aL (l2005 b it has been shown that the redshift 
distribution of flag-2/4 objects is different from that of the flag- 
0/1 objects. If this is also the case here, the above correction 
may not be accurate. The situation is less clear in the redshift 
range 3 < z < 4 since the CRFs of both flag-1 and flag-2 
objects are comparable to that of flag-1 objects at low redshift. 
In any case, this correction is probably valid only to the first 
order, and, in the following, we shall compare the LFs with and 
without this correction, which we call respectively "corrected" 
and "minimal" respectively. 

Flag-1, and to a lesser extent flag-2, objects may have been 
wrongly assigned a redshift outside of the range 3 < z < 4, 
and a full incompleteness correction should include a contribu- 
tion from these flag-1 and flag-2 objects erroneously classified 
as low-redshift galaxies. However, such correction is more dif- 
ficult to determine in a robust way, and we prefer not to take 
these objects into account, keeping in mind that the true galaxy 
counts are probably higher than those we calculated. 

3.5. Cosmic variance 

Cosmic variance, i.e. the presence of inhomogeneities due to 
the large-scale structure of the Universe, can play a role in the 
estimation of observables like surface densities and LFs, in par- 
ticular when one compares the results from different studies. 
The usual strategy to evaluate this effect is to use different fields 
and to analyze them separately (e.g., S99; ST06A). 

The VVDS Deep Field consists of a single field, but we can 
have a grip on the amplitude of cosmic variance by dividing the 
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Fig. 2. Surface density of VVDS high-redshift objects in sub- 
fields of different sizes compared to that of the full field (from 
left to right, 3x3, 2x2, 1x2 and 1x1). The statistical uncertainy 
on the full field is indicated. The grey boxes show the expected 
statistical uncertainties for the smaller fields. Number counts 
are corrected for the estimated fraction of correct objects. 

field into smaller fields. Fig. [2] show s the surface density of the 
high-redshift galaxies discussed in iLe Fevre et al. (2005a) in 
the VVDS Deep Field, together with the surface densities cal- 
culated on sub-areas obtained by cutting the full field in 2, 4 and 
9 smaller fields respectively. The purely statistical uncertainty 
on the full field can be extrapolated to the smaller fields by as- 
signing the galaxies at random to the different subfields. Fig.|2] 
shows that the dispersion in the smaller fields is consistent with 
that expected from pure counting statistics on the galaxies in 
the full field, and therefore that any effect of cosmic variance 
in the VVDS field on these scales is quite small compared to the 
statistical uncertainties. This does not preclude that structures 
on scales significantly larger than the full VVDS field may ex- 
ist; however, these structures are not predicted by cosmological 
simulations dSomerville et al.ll2004l) . In the following, cosmic 
variance is assumed to be negligible and is therefore left out. 



4. Results 

4.1. Ultraviolet luminosity function 

We calculate the luminosity functio n using ALF. A LF incorpo- 
rates the classical estimato r 1 /Vmax (ISchmidtll968h . the C"^ es- 
timator dZucca e t al. 1997), the maximum-likelihood estimator 
SWML dEfstathiou et al.„1988i) and the param etric maximum- 



1979I) . which uses the 



likelihood estimator ST Y dSandage et al 
distribution proposed bv lSchechteildl97d) 

0(M)dM = 0.4^*ln(10)(lO''-4(^*-^)"'^'exp(-10°-4(^'-^)dM 

In the discussion below, we shall focus on the standard pa- 
rameters obtained from STY: the characteristic absolute mag- 
nitude M*, the characteristic density (p* (in mag"' Mpc"^) and 
the slope a. 

To allow comparisons with previous studies (e.g., S99, 
ST06A), we calculate the luminosity functions in a ultravio- 
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Fig. 3. Luminosity functions for VVDS galaxies with 3 < z < 4 using the equal-weight method (a) and the photometric-rejection 
method (b). Top panels: 1700A luminosity function as a function of Mnoo- a has been fixed to -1.4. The heavy black and 
grey lines are the corrected and "minimal" LFs (see Sect. 13.4b . respectively. The empty symbols show the three non-parametric 
estimators for the corrected LF: 1 / Vmax (squares), (triangles) and SWML (circles). The black and grey squares are the 1 / Vmax 
estimators from ST06A at z ~ 3 and z ~ 4- respectively, and the thin lines are the corresponding STY LFs. The black stars are 
the non-parametric estimators of the z ~ 3 UV LF from IPoli et al. (2001). Bottom panel: Uncertainty contours of the 1700A 
luminosity function in the M*-(p* plane. The filled circles and the contours show the VVDS LFs; the minimal LF is shown in 
grey and the corrected LF in black. The heavy and thin solid lines are the 68% and 90% contours. Stars show the STY parameters 
obtained in other studies (see text); black, empty and grey stars are for z ~ 3, z ~ 3.5 and z ~ 4 respectively. The 68% uncertainty 
contours from ST06A are also shown. 



let pseudo-filter centered around 1700 A. The pseudo-filter we 
used is a 100 A- wide box centered on 1700 A, but the results 
do not depend on the details of the filter. The redshift range 
3 < z < 4 is particularly adapted to our observational set- 
up. At these redshifts, 1700 A falls very close to, or even in- 
side, the I filter used to select our galaxies. Absolute lumi- 
nosities, which are calculat ed by fitti ng theoretical templates 
from the PEGASE2 library jFioc & R occa-Volmerange 1997) 
to the full set of the nine photometric filters from the CFH12K 
and CFHTLS observations, will therefore need very small k- 
corrections. In addition, the redshift dist ribution of the objects 
in the VVDS falls quickly below z ~ 3 dLe Fevre et alJl2005bl) 
as a result of the so-called "redshift desert", produced by the 
combination of the observed wavelength range with VIMOS 
and the paucity of spectral lines and features in the rest-frame 
UV 1800-3500 A galaxy spectra. The spectroscopic nature of 
the VVDS also Hmits the depth of the survey above z ~ 4. 



In the redshift range 3 < z < 4 the VVDS reaches abso- 
lute magnitudes Mnoo ^ -21.5, which makes the a param- 
eter completely unconstrained, as this limiting magnitude is 
brighter than M* found in previous studies of the UV LF. Other 
attempts to derive the high-redshift UV LF found rather con- 
flicting values of a, ranging from an assumed value of - 1 .6 in 
S99 to ~ -1.4 in ST06A, a nd even to ~ - 1.07 in the FORS 
Deep Field study (FDF; Gab asch et alj|2004 . Instead of trying 
to determine the value of a, we fixed it to -1.6, -1.4 and -1.1 
to evaluate its effect. As most recent studies found values com- 
patible with -1 .4 (see Sect. 15. lb . we shall consider this value of 
a more closely. The resulting equal-weight LFs with a = - 1 .4 
are presented in Fig. [3^, and the photometric -rejection LFs 
in Fig. [3J3. The non-parametric estimators are quite consistent 
with each other even in the lowest-luminosity bin, which in- 
dicat es that there is little bias in the population under study 
(Ilber t et al.l 120041) . Uncertainty contours in the M*-(f)* plane 
for a - - 1 .4 are plotted in the bottom panels of Figs[3^ and[3}3. 
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Table 1. Parameters of the UV luminosity functions estimated 
from the VVDS for different fixed values of a. "EW" indicates 
the equal-weight LPs and "PR" the photometric-rejection ones. 
(p* is in units of 10"^ mag^' Mpc^^ The 68% and 90% (in 
parenthesis) confidence levels are also provided. 



Method (/>' M' a 

1.6 
1.4 
1.1 

1.6 
1.4 
1.1 

1.6 
1.4 
1.1 

1.6 
1.4 
1.1 



Table 2. UV luminosity density -Cnoo(M < -18.5) integrated 
up to Mi7()o = -18.5 and total £J°^^^ estimated from the VVDS 
for different fixed values of a. "EW" indicates the equal-weight 
LFs and "PR" the photometric -rejection ones. The 68% and 
90% (in parenthesis) confidence levels are also provided. Units 
are 10"^ WMpc^^ 



Method i;,7oo(M < -18.5) £]f^ a 



"EW" minimal 


13.88!|^^ 
11.78!;;^^ 

y.KJJ Qgg 


■+2.76^ 
-2.37, 
f +2.251 
1-1.90) 
f +1.841 
1,-1.47; 


19 20+"-'' 
13 91+^-'^ 


^+4.47^1 
1-3.791 
f +3.001 
1-248 J 
I +2.061 
1,-1.64; 


-1.6 
-1.4 
-1.1 


"EW" corrected 


17 54+^ " 
14 84+1 68 

12.02!];! 


^+3.94^1 
1-3.38 J 
f +3.031 
1-2.67; 
1+2.331 
1-2.00; 


24 75+3-52 
12.90!|:« 


f+6.45"! 
1-5.33) 
1 +4.09 1 
1-3.51) 
1 +2.62 1 
1-2.25; 


-1.6 
-1.4 
-1.1 


"PR" minimal 


22.15!|« 
18.45^3 g 
14.64!|^2 


^+9.891 
1-6.88; 
f +8.381 
1-5.61; 
("+5.60 
1-4.27; 


33.50!^;"° 

i6-ii!3:« 


r+17.5^ 
1-11.9) 
f + 11.7' 
l-V.78; 
("+6.60 
1-4.95; 


-1.6 
-1.4 
-1.1 


"PR" corrected 


29.20!^-^ 
24.31^1^2 

18.98!3;^8 


+ 13.2"! 
1-9.85; 
('+11.0~1 
1-7.93, 
f +7.641 
1-5.51; 


45.07!;^;^ 
30.96!?:^^ 
21.04!^:« 


^+24.0 
1-17.1; 
r+15.3^ 
1-11.1, 
1+9.181 
1-6.44; 


-1.6 
-1.4 
-1.1 



Q 26+0-'" (+"131 -22 22 

"•^"-0.07 1-0.10/ zz,.z.z 

"EW" minimal 0.32+° j;^ H fA -22.07 

0-38t««(t!^) -21.89: 

o.38tlS(!S!:?§) -22.08: 

"EW" corrected 0.45+j; |[ (^jS J^j -21.95 

0.53tl|(t??) -21.78: 

0.75t^H-o:I^) -21-71 

"PR" minimal 0.85!;;:3^ (^Sj;^') -21.59: 



n Q4+0.33 /+0.70) aa 

"■^^-0.34 (-0.46,1 ^ ^ 

1 14+0.49 ('+1.06 1 _9i cq 
^•^^-0.50 1-0.66) 

"PR" corrected 1-24!;; *|5 -21.49: 

1 Q4+0.48 /+1.02\ -ir 

-0.49v-0.66y ^i.JJ 



+0.16 


f+0.26' 


-0.17 


1-0.28 


+0.15 


1+0.23 


-0.15 


1-0.25, 


+0.13 


+0.21 


-0.13 


-0.22, 


+0.15 


f+0.24' 


-015 


l-0.26i 


+0.14 


f+0.22\ 


-0.14 


1,-0.23) 


+0.12 


f+0.20' 


-0.12 


1-0-21) 


+0.22 


f+0.35\ 


-0.22 


^-0.38) 


+0.20 


(+0.33\ 


-0.21 


l-0.35i 


+0.18 


f+0.29 


-0.18 


^-0.30) 


+0.20 


(+032' 


-0.20 


1-0.35 


+0.19 


f+0.30 


-0.19 


1-0.32 


+0.17 


f+0.27 


-0.17 


1,-0.29, 



The Schechter parameters from the STY estimator are listed in 
Table □ 

The choice of the method (equal-weight vs photometric- 
rejection) has a significant impact on M*. Other parameters be- 
ing equal, its value is about 0.5 magnitude brighter in the equal- 
weight LFs than in the photometric-rejection ones. This is con- 
sistent with the rejection by the latter method of low-redshift 
galaxies that were mistakenly classified as high-redshift galax- 
ies. As the density drops rapidly below M*, even a small num- 
ber of these spurious galaxies may alter significantly the esti- 
mate of M*. Therefore we consider the photometric -rejection 
LFs to be more robust and adopt the results from this method 
as our best estimates. 

The correction for flag-0 objects increases (p* by about 40- 
50%, but also increases M* by about 0.1 mag, both effects be- 
ing expected, since the correction gives more weight to the 
faintest galaxies in our sample. This correction is quite con- 
servative, because an additional contribution from flag-1 and 
flag-2 objects erroneously classified as low-redshift galaxies is 
to be expected in addition to the flag-0 objects (see Sect. 13.4b . 
However, this additional correction, which would also increase 
(p*, is probably small and less straightforward to include in a 
robust way, and we did not attempt to include it. 

The error contours in Figs [3^ and[3}3 show that, for a fixed 
slope a - -1.4, M* and (p* are quite strongly correlated, a 
bright M* implying a low (p* and vice versa. The same ef- 
fect can be seen by varying a, with a change on M* of about 
+0.15 mag between a - -1.1 and a = -1.6. 

In conclusion, when fixing the slope to a = - 1 .4, the best 
estimate of the UV 1700 A LF in the range 3 < z < 4 using 
the VVDS spectroscopic survey has the following STY param 
eters: cp* = 1.24+j;-^g 10"^^ mag"' Mpc^^ and M* 
(68% confidence levels). 



4.2. Ultraviolet luminosity density 

The luminosity density (LD) at 1700 A -Cnoo, which is the to- 
tal power per volume unit contained in 1700 A rest- wavelength 
photons emitted by galaxies, can be calculated by summing, 
possibly up to a limiting magnitude Mum, the contributions of 
the individual galaxies weighted by their number density given 
by the LF: 



00 



-1700 



(M) (p(M)dM, 



(1) 



21 49+"- 19 



where Li7oo(Mi7oo) = IQ-^-'^iMnw-M.i) jj^-i luminos- 
ity at 1700 A of a galaxy with an absolute magnitude Mnoo- 
The luminosity density X is expressed in W Hz ' Mpc"^; the 
total LD is the limit 00. 

The major difficulty in the calculation of the luminosity 
density by integrating the LF is the effect of the a parame- 
ter of the Schechter distribution. Indeed, for steep values of 
a, the contribution of the faintest galaxies, where the LF is 
the least constrained, may dominate that of the much less nu- 
merous bright galaxies. The usual work-around is to adopt a 
limiting absolute magnitude Mum, although its choice is non- 
physical and rather arbitrary. Even when such choice is war- 
ranted, any uncertainty in a may result in significant changes 
of the LD. We address this problem by calculating the LDs for 
different fixed values of a's. We follow ISawicki & Thompson 
(2006b), who calculated, using the LFs in ST06A, the LD up to 
^1700 = -18.5, as well as the total LD. The results are shown 
in Figs|4^ and|4j). 

Expectedly, the LDs calculated from the uncorrected LFs 
are systematically lower than those from the corrected ones, 
the amplitude of this effect being between 20 and 50%. When 
integrating up to Mnoo = -18.5, the effect of changing a 
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Fig. 4. Luminosity density at 1700 A integrated up to M = -18.5 (a) and total (b) as a function of the assumed value of a, based 
on the photometric -rejection LFs. Circles show the LDs including the flag-0 correction, while squares do not. Stars show results 
from other studies (see text). Black, empty and grey stars are for z ~ 3, z ~ 3.5 and z ~ 4 respectively. Symbols at a = -1.6 have 
been sUghtly shifted for better readability. The different lines connect related symbols for easier identification only. 



is moderate, as the three LDs are consistent with a value 
of Xi70()(M < -18.5) ^ 23 10'^ W Hz"' Mpc^^ for the 
photometric-rejection LFs. This value is about 2.5 times larger 
than -C\ioo{M < -21.0), which is the limit of integration used 
in S99 to avoid extrapolating the LFs. The density of UV pho- 
tons from galaxies in the range -21.0 < Mnoo < -18.5 
is therefore comparable to that of the brightest galaxies with 
MnoQ < -21.0. When integrating the full LD, the effect of a 
becomes quite serious, as the LFs differ by about a factor 2 be- 
tween the a - -1.1 and a - -1.6 cases. Except for the very 
steep a = -1.6 case, £noo(M < -18.5) is > 80% of 
meaning that the bulk of the UV photons comes from galaxies 
brighter than M* + 3, i.e. 15 times fainter than the characteristic 
luminosity. 

5. Discussion 

5.1. Comparison of the LF with previous studies 

This study of the ultraviolet LF at 3 < z < 4 is the first one 
entirely based on a spectroscopic sample. It is therefore impor- 
tant to compare the results in detail with prev ious works using 
LBG color selecti on rS99:IOuchi et al.ll2004 ST06A) o r pho- 
tometric redshifts ("Poli et al."200lVGabasch et al."2004. 

With the exception of the KDF sample (ST06A) and the 
HDF sample (IPoli et al.ll200lh . some work was needed before 
we could compare our results with the previous ones, because 
of differences in the calculation of the LFs. S99 give their STY 
parameters for a cosmology with Hq - 50 km s ' Mpc ' and 
Qfn - 1; transforming the STY parameters of S99 into the 
concordance cosmology, we get M* = -21.07, (p* - 1.4 lO"-' 
3andM* = -21.12, £^ 



mag ' Mpc ^ at z 

Mpc"^ at z ~ 4. In the case of the FDF (Gaba sch et al 



1.07 10-^ mag-' 



2004) 



a pseudo-filter at 1500 A has been used; we therefore corrected 
their value of M* by -0.1, which corresponds to the difference 
in M* we obtain if we calculate the VVDS LFs at 1500 A in- 
stead of 1700 A. In the case of the SDS dOuchi et al.,,2004.) . 



which found a much steeper a than all other groups, we chose 
their STY parameters where a has been fixed to -1 .6. 

The above surveys are photometric (howe ver, the S99 



sample has been spectroscopically confirmed in ISteidel et al 
(1996)), which allows them to reach deeper magnitudes than 
the VVDS, which is a pure magnitude-limited spectroscopic 
survey. The shallower depth of the VVDS doesn't allow us 
to constrain the slope a independently. We therefore fixed the 
value of a to different values found in the literature, and fi- 
nally adopted a = -1.4, compatible with the latest and deep- 
est study of the LF of LBGs (ST06A), as we ll as with studies 



based on photomet ric redshifts in the HDF (IPoli et al .112001 



Arnouts et al.ll2005h . It cannot be taken for granted, however. 



that the slope we would have found with a much deeper spec- 
troscopic survey would be compatible with those obtained with 
Lyman-break or photometric -redshift selections, for instance if 
these selections miss a significant fraction of the high-redshift 
galaxies (see Sect. 15. 2b . In any case, the discussion of the value 
of a is beyond the scope of this paper 

The value we obtain for (p* for the PR LF corrected 
for incompleteness is compatible with those found in stud- 
ies of the LBG pop ulation (S99; ST06A), except for the SDS 



(lOuchi et al.l 120041) . whose LF at z ~ 4 is extreme compared 



to all other studies (see Fig. |3] bottom panels). Photometric- 
redshift studies of the HDF and FDF find however slightly 
larger values of (p* which are marginally compatible with ours 
at the 90% level. We obtained our result by applying only one 
very simple correction for incompleteness that considers only 
flag-0 objects. A correction for missed 3 < z < 4 objects which 
ended up with flags 1 and 2 in a different redshift range would 
increase 0*, but we conclude that we do not find any evidence 
in the VVDS high-redshift LF for a value of 0* in significant 
disagreement with values found in previous studies. 

With a value of M* - -21.49, the characteristic luminos- 
ity is found to be ~ 0.5 mag brighter than that found in other 
studies, which are close to M* = -21.0, or even fainter (SDS). 
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Our M* estimate differs from these results at the 2.5cr level at 
least. We note that M* may even be brighter if the photomet- 
ric rejection erroneously discards some very luminous galaxies 
with correct redshifts (see the values of M* found for the EW 
LPs in Table [TJ. 

As shown in the bottom panels of Fig. [3] the M*-4>* error 
contours are strongly correlated. If M* has been underestimated 
because of statistical fluctuations, (p* would also have been un- 
derestimated, pushing it to values quite higher than those found 
in S99 and ST06A at the level of the HDF (" Poh et al.ll200ll) 



and the PDF dGabasch et alj|2004.) : however, even pushing the 
statistics to its limit, it remains difficult to reconcile the charac- 
teristic absolute magnitude M* in the VVDS high-redshift LP 
with those obtained from other LPs. We discuss this further in 
the next section. 




5.2. A new population of non-LBG high-redshift 
galaxies? 

Assuming a - -1.4, M* is about -2L5, which is 60% brighter 
brighter than all values of M* obtained so far at this redshift 
and this wavelength. This quite different M* needs to be under- 
stood. 

While Poisson statistics is included in our uncertainties, 
we checked that cosmic variance on scales commensurable 
with the size of the VVDS Deep Pield is probably negli- 
gible compared to our statistical uncertainties. Structures on 
much larger scales may exist, but are not predicted by theory 
dSomerville et al.l 120041) . Cosmic variance affecting the fields 
used by other authors cannot probably explain the differences 
in the LPs. The KDP survey of ST06A covers about a tenth 
of the area of the VVDS Deep Pield (169arcmin^ compared 
to 1 720arcmin^) split into five fields. The effect of the cosmic 
variance on (p* was found to be of the order of 50%. 

The reasonably good agreement between most previous 
studies is expected, since they are based on similar observ- 
ing strategies, and therefore may be subject to similar bi- 
ases. S99 and ST06A used the same LBG-selection strat- 
egy to identify Lyman-break galaxies, cle arly different from 
our pure magnitude-limited selection; SDS ( Ouchi et alj|2004 ) 
also uses an LBG selection, b ut with different filters. Studies 
based on photometric redshifts ( Poli et al1l2001 ; Gabasch et al 



-0.5 



0,5 



1,5 



2004l:lArnouts et alj|2005b heavily rely on similar assumptions 
on galaxy emission and intervening absorption. It is possi- 
ble that these assumptions do not apply to the entire high- 
redshift galaxy population. Furthermore, the presence of a 
strong break in the bluest band remains the most important cri- 
terion in identifying a galaxy as a high-redshift object, making 
the photometric -redshift population little different from that of 
LBGs in practice. It must be pointed out that straight galaxy 
counts are much larger at all magnitude s in the VVDS than, for 
instance, in S99 ( he Pevre et al.ll2005a ^. and that the relatively 
good agreement of (p* with LBG-based studies is due to the fact 
that the latter include substantial incompleteness corrections 
based on complex simulations that are completely avoided in 
the VVDS. 



Fig. 5. (u* - g') vs (g' - r') color-color diagram for the VVDS 
galaxies with redshifts in the range 3 < z < 4. Upper limits on 
u* are indicated by a vertical arrow. Black squares are objects 
with secure redshifts (flags 3 and 4). Empty squares and crosses 
are objects with less secure redshifts (flags 1 and 2). A square 
indicates an object kept in the calculation of the PR LP. The 
solid line shows the contour of the Lyman-break selection box. 
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Fig. 6. Composite spectrum of all 172 3 < z < 4 galaxy candi- 
dates which fall outside of the LBG selection box. 



The VVDS is essentially free of any assumption on the 
presence of a break in the continuum; the photometric-rejection 
LPs reject galaxies with very improbable colors, which is not 
at all the same as selecting galaxies based on their colors. The 
bright M* found in our survey may therefore be an indication 
that there exists a population of high-redshift galaxies with "un- 
expected" colors. Pig. |5] shows the (m* - g') vs (g' - r') color- 
color diagram for our 3 < z < 4 galaxies. The Lyman-break 
selection box has been calculated on the basis of the criteria 
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of ISteidel et aD d 19961) and adapted to the CFHTLS filters, so 
that it selects efficiently galaxies with redshifts 3 < z < 4; we 
have adopted a quite conservative bottom limit for the box, so 
that it should contain the majority of the LBGs, at the price 
of some contamination by low-redshift galaxies. Among the 
twelve secure-redshift galaxies, only four of them fall inside 
the selection box; three more objects, having only upper limits 
on u*, may also be inside the box. Five galaxies are however 
clearly outside the box. Incidentally, Fig. |5] shows that most 
of the objects rejected for the calculation of the PR LF have 
g' - r' < 0.4, a region where numerous low-redshift galax- 
ies are expected, which confirms the validity of the method. 
Fig. |6] shows the composite spectrum of the 172 galaxies in 
our 3 < z < 4 sample whose colors place them outside of 
the LBG selection box. The composite spectrum presents very 
clearly all the expected features of galaxies at z ~ 3-4, both 
in the spectral lines and in the continuum, making the case for 
the existence of numerous high-redshift galaxies outside of the 
selection box quite convincing. One would expect that the LF 
of VVDS galaxies inside the LBG selection box would match 
those from S99 and ST06A; however, the calculation of such 
LFs is quite complex, because it requires the understanding and 
the modelling of all effects that may result in a failure in the 
magnitude determination in any of the three filters, an accurate 
estimate of the contamination by low-redshift galaxies, as well 
as the "natural" crossing of the box boundaries because of the 
statistical uncertainties in the magnitudes. Such work is there- 
fore outside of the scope of this paper. 




Fig. 7. Li7O()(M)0(M)/X'[°Qg as a function of M for our best 
estimate of the LD (black line) calculated from the PR LF, 
compared to the LD at z ~ 3 from the KDF (heavy grey line) 
and at z ~ 4 from the SDS (thin grey line). The black thin 
fines show the 68% uncertainties on Li7oo(M) 0(M)/X|°Qg. The 
dashed vertical line indicates M - -21.0. 



One possible reason for the presence of galaxies out- 
side of the LBG selection would be that the intergalactic ab- 
sorption has been overestimated. Indeed, by recalculating the 
transmission t h rough the intergalactic medium first derived by 
Madaul (Il995h . iMeiksinI (12006) found differences in colors of 
about 0.5-1 mag for bands containing rest-frame Lyman a and 
shorter wavelengths. Should the bottom of the selection box 
be lowered by about 0.5 mag, as suggested by figure 2b in 
Meiksin (200a), it would become hard, if not impossible, to 
use LBG selections without spectroscopic follow-ups to build 
complete samples of high-redshift galaxies, as contamination 
by low-redshift galaxies will be important in the lower part of 
the box. Particularly strong absorption in the G filter, which is 
affected by the presence of the Lyman a forest, might also con- 
tribute to the presence of high-redshift galaxies outside of the 
LBG selection box. 

At redshift z ~ 2, it is known that BM/BX-selected galaxies 
(whi ch is an extension at low er redshift of the LBG color cri- 
teria; lAdelberger et al.ll2004l) do not form the complete popu- 
lation of star-forming galaxies. Combinin g different photomet- 
ric selection criteria, Reddv et al. ('2005") found that the star- 
forming BzK galaxies (Daddi et al.. 2004) do not all satisfy 
the LBG color criteria, the overlap being only about 70-80%. 
The BzK color criterion is more efficient at selecting IR-bright 
galaxies, and these galaxies might have redder G - H colors 
than the BM/BX-selected ones. The VVDS might therefore ob- 
serve the equivalent of the BzK star-forming galaxies at red- 
shifts 3 < z < 4, in addition to the "normal" LBG galaxies. 



5.3. The ultraviolet luminosity density In the VVDS 

We compare the UV LD obtained from the VVDS LF at red- 
shifts 3 < z < 4 with those obtained in the studies discussed 
in Sect. 15.11 Most of these studies quote the values of the LD, 
but, when needed, we calculated the LDs ourselves through the 
integration of the Schechter distributions they obtained. 

In the calculation of the LD, the choice of a plays a very 
important role, making comparisons between LFs with differ- 
ent a quite meaningless. This effect is quite visible in Figs|4^ 
and (especially) lit. Because of the large error bars inherent in 
the determination of the LDs, Iinoo{M < -18.5), i.e. the LD 
integrated up to Mnoo - -18.5, and X}°Jfo. the total LD, are 
mostly compatible with the resu lts from previous works; only 
the LD at z ~ 4 from the KDF dSawicki & Thompsonll2b06bl) 
being markedly lower. 

While the total LD from the VVDS is not very different 
from the already published values, the significantly brighter M* 
we have found modifies strongly the picture of the distribution 
of the LD as a function of the galaxy luminosity. The VVDS 
LD integrated up to Mnoo - - 21.0 is indeed more than twic e 
that from the KDF at z ~ 3 dSawicki & Thompsonir2006bl) ; 
this is akeady seen in Fig. |3j>, where our non-parametric LF 
is always at least twice as large as that of ST06A. Galaxies 
brighter than Mnm - -21.0 account for a third of the total LD 
in the VVDS, while this figure is only 15% in the KDF The 



same appfies if we compare with other studie s (S99; lPoli et al 



20011: lOuchi et al.ll2004t iGabasch et al.ll2004 . Adopting com- 



patible values for a, the ratio -Ci70o(M < -21.0)/X|™ is 
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typically twice as large in the VVDS compared to the other 
studies. We conclude that the VVDS shows evidence of a 
more important contribution of the most UV-bright galaxies 
to the total LD. This is illustrated on Fig. |7] which shows 
i'i70o(A/)</'(-^)/-C|°oo, i.e. the fraction of the luminosity den- 
sity at magnitude M to the total LD (see Eq. ([1])), as a function 
of absolute magnitude M for our study and for the KDF study. 
Comparisons with other studies qualitatively agree, the ampli- 
tude of the effect being the strongest for the SDS. 

6. Conclusion 

The UV 1700 A luminosity function in the VVDS for galax- 
ies with redshifts 3 < z < 4 differs significantly from other 
estimates of the UV LF in the value of the characteristic lumi- 
nosity, above which the number density of galaxies decreases 
exponentially; this characteristic luminosity, corresponding to 
an absolute magnitude Mnoo ~ -21.5, is found to be about 
60% larger than previously reported. On the other hand, the 
galaxy number density at lower luminosities is mostly compat- 
ible with previous studies. Our result has been obtained with 
minimal corrections for incompleteness, and it is possible that 
the remaining incompleteness would significantly increase the 
total luminosity density. 

The difference in the LF seems to be due to the fact that 
Lyman-break selection techniques and photometric -redshift 
studies are not able to retrieve the full population of high- 
redshift galaxies. This may be due to a combination of differ- 
ent effects, like, for example, photometric problems that are not 
well modelled, intergalactic absorption that pushes the bottom 
of the Lyman-break selection box down in the color-color di- 
agram, or a population of high-redshift galaxies whose colors 
do not correspond to the expectations of Lyman-break galax- 
ies or photometric redshifts. This would make color-selection 
techniques unable to build samples of high-redshift galaxies 
representative of the complete galaxy population. The popu- 
lation of high-redshift luminous galaxies may include a mix 
of LBGs, the high-redshift equivalent of the BzK star-forming 
galaxy population, and possibly other types of galaxies. 

The determination of the luminosity density requires a pre- 
cise determination of the slope a of the luminosity function 
at faint magnitudes. Even though some very deep photomet- 
ric studies like the KDF (ST06A) have put strong constraints 
on a, it is not clear whether the same a would be found if the 
full population of high-redshift galaxies could be included. In 
absence of color criteria able to recover the entire population, 
the determination of a will have to await spectroscopic surveys 
similar to the VVDS but reaching significantly fainter magni- 
tudes. 

While we find that the total luminosity density in the VVDS 
is compatible with previous estimates, the distribution of the 
LD among galaxies of different luminosity is quite different, 

with galaxies brighter than Mnoo 21.0 accounting for about 

a third of the total LD, a fraction twice as large as previously 
estimated. Therefore, the VVDS paints a quite different picture 
of the role of the most actively star-forming galaxies in the his- 
tory of star formation. 
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